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Abstract 

In their simplest formulations, classical dynamics is the study of 
Hamiltonian flows and quantum mechanics that of propagators. Both 
are linked, and emerge from the datum of a single classical concept, the 
Hamiltonian function. We study and emphasize the analogies between 
Hamiltonian flows and quantum propagators; this allows us to verify 
G. Mackey’s observation that quantum mechanics (in its Weyl formu¬ 
lation) is a refinement of Hamiltonian mechanics. We discuss in detail 
the metaplectic representation, which very explicitly shows the close 
relationship between classical mechanics and quantum mechanics, the 
latter emerging from the first by lifting Hamiltonian flows to the dou¬ 
ble covering of the symplectic group. We also give explicit formulas for 
the factorization of Hamiltonian flows into simpler flows, and prove a 
quantum counterpart of these results. 
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1 Introduction 


Hamiltonian and quantum mechanics share a common background: they 
are both based on the datum of a function H (the Hamiltonian) which is 
a generalization of the notion of energy. While this function is taken liter¬ 
ally in Hamiltonian mechanics as defining a vector field Xh on phase space 
which gives rise to the “Hamiltonian flow” in quantum mechanics one 

associates with H a differential, or pseudo-differential, operator generating 
a group of unitary transforms (U^)t satisfying the Schrodinger equation. 
When the Hamiltonian function is a quadratic form in the phase space vari¬ 
ables, both and (U f H )t are easily linked: in this case the f t H are just 

linear canonical transformations, and the Uf 1 are (up to an easily determined 
phase factor) the metaplectic operators obtained by lifting the flow 
to the metaplectic group. We will, among other things, extend this simi¬ 
larity to arbitrary Hamiltonian flows and propagators, thus producing an 
extension of the metaplectic representation. One of the aims of the present 
paper is actually to highlight the similarities between classical mechanics (in 
its Hamiltonian formulation), and quantum mechanics: the latter is seen as 
emerging from the first. In fact, all which lacks is a physical motivation for 
the introduction of Planck’s constant h. We will not discuss this very im¬ 
portant - and indeed difficult! - problem here; for all practical purposes we 
will view Planck’s constant as a scaling factor, and we will not discuss the 
physical content of quantum mechanics (in the Schrodinger picture, which 
is the only one considered here): there is a vast literature on the subject, 
and it seems that no real consensus has yet been reached. 

A point of terminology: we will indifferently use the words “canonical 
transformation” or “symplectomorphism”. While in some texts both no¬ 
tions are not quite equivalent, “canonical transformation” sometimes hav¬ 
ing meaning any diffeomorphism preserving the form of Hamilton’s equa¬ 
tions, we will adhere to the stricter definition following which a canonical 
transformation is a diffeomorphism of phase space whose Jacobian matrix 
is symplectic. 

Here are some highlights of this paper: 

• In Section [2] we study Hamiltonian flows from the symplectic point of 
view, and prove factorization formulas. For instance we show in Propo¬ 
sition [7] that if a Hamiltonian H is split into a sum of two independent 
Hamiltonians Hq + H\ then the flow (ff H )t determined by H can be 
written as the product f f H ° //^ 1 where H[(z,t) = H\ ( f t H ° (z). t). These 
constructions allow us to define the group Ham(n) of all Hamiltonian 
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symplectomorphisms associated with compactly supported Hamilto¬ 
nians. We also discuss a few special properties of Hamiltonian flows 
(and canonical transformations) such as the theorems of Katok and 
Gromov. 

• In Section [3] we give a rigorous review of the metaplectic group and 
of its inhomogeneous extension. This section is in a sense pivotal, 
because it is the gate to quantum mechanics, but a gate only using 
classical concepts (the path lifting property for covering groups). It 
shows, in a sense, that quantum mechanics emerges from classical 
(Hamiltonian) mechanics by “passing to the covering group” (there is 
a certain similarity with the passage from the rotations group to the 
spin group). We discuss the notion of Feichtinger algebra which is the 
most natural domain for the metaplectic group and its inhomogeneous 
extension. We finally shortly review the theory of the Weyl symbol of 
a metaplectic operator, following previous work of ours. 

• In Section 01 after having discussed quantization in general, we de¬ 
fine the notion of quantum isotopy: a quantum isotopy is a continu¬ 
ous path (Ut)t of unitary operators satisfying a certain differentiabil¬ 
ity condition which allows us to define a canonical self-adjoint opera¬ 
tor by the formula H = ih (4f Ut) Uj~ . This operator coincides with 
the infinitesimal generator obtained by Stone’s theorem when ( Ut)t is 
a strongly continuous one-parameter group. We thereafter prove in 
Proposition 03 a quantum analogue of the factorization result of [TJ we 
have U?° +Hl = U t H ° U^ 1 oSt if Hq is quadratic. 

Notation We will identify T*W l = M n x (R n )* with M 2n and denote by 
a the standard symplectic 2-form J2j=i^Pj A dxj, i.e. a(z,z') = Jz ■ z' 

where z = (x,p), z' = (. x' ,p') and J = ^ ^ ^ . The scalar product of two 

vectors u, v £ M m is written u ■ v, if A is a symmetric m x m matrix we will 
often write Au ■ u = Au 2 . We will denote by h a positive parameter which is 
identified in physics with Planck’s constant h divided by 27r. The Schwartz 
space of complex functions on W n decreasing at infinity, together with their 
derivatives, faster than the inverse of any polynomial is denoted by 5(M m ); 
its dual S'(W n ) is the space of tempered distributions. The scalar product 
on L 2 (M m ) is written (-|-). 
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2 Symplectic and Hamiltonian Isotopies 


The symplectic group Sp(n) is the group of all (linear) automorphisms of 
the symplectic space (R 2n ,<r): S £ Sp(n) if and only if S £ GL(2n,R) and 
c t(Sz,Sz') = a(z,z') for all (z, z') £ R 2n x R 2n (for a detailed exposition, 
see |28|). A diffeomorphism / of R 2n is a symplectomorphism if its Jacobian 
matrix at every point is symplectic: Df(z ) £ Sp(n) for every 2 : £ R 2n . 
Equivalently, = a, identifying a with the 2-form Y^j=i^Pj A dxj. 

The symplectomorphisms of (R 2n ,<r) form a group Symp(n). 

2.1 Hamiltonian vector fields and flows 

2.1.1 Hamilton’s equations 

Let H £ C' j (R 2ti x R), j > 2, be a real-valued function; we will call H a 
Hamiltonian function. The associated Hamilton equations with initial data 
z’ at time t' are 


x = d p H(x,p,t) , p = -d x H(x,p,t) (1) 

or, using the collective notation z = (x,p), 

z{t) = X H (z(t),t) , z(t!) = z ' (2) 

where Xh = Jd z H is the Hamilton vector field (strictly speaking it is not 
a true vector field when H depends on t). Existence and uniqueness results 
for Hamilton’s equations abound in the literature (see e.g. EI21EZ1 and 
the multiple references therein). The study of these properties are actually 
mostly a branch of the theory (local, and global) of dynamical systems 
(= systems of ordinary differential equations) where the notion of vector 
fields and their integral curves play the primordial role. The main result we 
will (implicitly) use is the following local classical existence and uniqueness 
property: Let X : R 2n —> R 2n be a vector field (Hamiltonian, or not) of 
class C J , j > 1 (it is hence, in particular, locally Lipschitz continuous). For 
every zo £ R 2n the system i = X(z) there exists an open ball B 2n (zo,r) C 
R 2n and an e > 0 such that this system has a unique solution 1 1 —> z(t) with 
z(0) = zo which is C J and defined for every t £ I e = [— s,e]. Moreover, by 
uniqueness, every such solution can be continued onto a so-called maximal 
interval It = [—T, T], The same result holds for time-dependent vector fields 
Xt , which are families of “true” vector fields depending in a C 1 fashion on 
time t. Transposed to the case of Hamiltonian systems, this means that it 
is sufficient, to have a local existence and uniqueness statement to assume 
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that H £ C- 7 (M 2n x M), j > 2 as we did. In practice, we will assume that 
the maximal interval It = [—T, T] is the same for every initial point zo, and 
whenever this simplifies statements, that T = +oo (which is the case if the 
Hamiltonian function is constant outside some compact set in M 2n , or, more 
generally, if the vector field Xjj is complete; see Abraham and Marsden [1] 
for details). 

We have the following transformation formula for Hamiltonian vector 
fields: if g £ Symp(n) then 

X H o g (z) = [Dg(z)]-'X H (g(z)). (3) 

This formula can be written more concisely as X gt h = g*Xu where g * is 
the operation of “pushing forward”. 

2.1.2 Hamiltonian flows 

Assuming existence and uniqueness of the solution for every choice of ( z',t ') 
the time-dependent flow (//(,) is the family of mappings M 2 " —> M 2ri which 
associates to every initial z' the value z(t) = f^ t ,(z') of the solution of 
We will write f t H = f t H () and we have 

= ftWfl)- 1 = ft H (f?)- 1 (4) 

and the f t H t , satisfy the groupoid property 

= ft H t» , ft?t = (5) 

for all t, t r and t". It follows from the first formula ([5]) that we have (//^,) _1 = 
f H 

h',t- 

The mappings f f H (and hence the mappings fll,) are symplectomor- 
phisrns, that is {ft 1 )*a = a. To see this, it suffices to notice that ix H & = dH 
and hence, in view of Cartan’s formula for the Lie derivative [1], 

Lx h ct = d(ix H cr) + ix H da = 0 

which implies that 

J t (ft H )*v = ft H YLx H CT = 0 

and hence {f^)*a = {f^)*cr = a (for a conceptually simpler, but longer, 
proof using Jacobian matrices, see [28]). 
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It follows from the transformation formula © for Hamiltonian vector 
fields that the flows (f t H )t and (f^° 9 )t are conjugate: 

ft ° 9 = <f l ft9\ (6) 

this is the property of “covariance of Hamilton’s equations under canonical 
transformations” familiar from physics (see Arnol’d [3], Hofer and Zehnder 
[55]). 

2.1.3 Special features of Hamiltonian flows 

Hamiltonian flows are volume preserving (because the Jacobian determinant 
of a canonical transformation is always equal to one). They however enjoy 
several unusual properties which makes them very different from arbitrary 
volume-preserving diffeomorphisms. For instance, Gromov proved in 1985 
that no Hamiltonian flow (or more generally, no canonical transformation) 
can embed a phase space ball B 2n (R) with radius R inside a cylinder Z 2n {r) 
based on a plane Xj,pj of conjugate coordinates if its radius r is smaller 
than R. This unexpected property, sometimes dubbed the “principle of the 
symplectic camel”, has led to various developments in symplectic topology. 
The principle of the symplectic camel has the following consequence: let us 
project orthogonally a ball with radius R on a plane of conjugate coordinates 
XjiPj, we then obtain a circle with area nR 2 . If we now deform the ball 
using a Hamiltonian flow, it will deform, but the area of the orthogonal 
projection of this deformed ball will never decrease below its initial value 
irR 2 . This result is a classical version of the uncertainty principle as we 
have explained in de Gosson [28, [31] [32] and de Gosson and Luef [42]. We 
mention that Gromov’s theorem in addition allows to define a new class 
of symplectic invariants, called symplectic capacities , which seem to play 
an important role in accelerator physics (see Dragt [14] and Erdelyi HZ])- 
The properties above are of a topological nature, and show that general 
volume preserving mapping cannot be approximated in the C° topology by 
canonical transformations. On the other hand, Katok m has showed that 
given two subsets fi and ft' with same volume, then for every e > 0 there 
exists a canonical transformation / such that Vol(/(H) \ O') < e. Thus, 
an arbitrarily large part of H can be symplectically embedded inside IT. A 
third property (which is however shared by other groups of diffeomorphisms) 
is IV-transitivity. Let us introduce the following terminology: given a group 
G acting on a set M we say that this action is IV-transitive if given two 
arbitrary sets {aq,..., xn} and {yi, ...,yi\r} of points of M there exists g € G 
such that g{xi) = yi for all i = 1,2, ..,1V. For instance, it is well-known 
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that Diff(M m ) (the group of diffeomorphisms of M m ) acts TV-transitively on 
for every TV. Now, a deep theorem of Boothby (9j (also see mm) 
says that the action of Ham(n) on M 2 ' 1 is TV-transitive for every integer TV. 
In fact, given two arbitrary sets {z\, ..., zn} and {z[, ..., z' N } of TV distinct 
points in M 2n , there exists a Hamiltonian function H such that z) = f^(zj) 
for every j = 1,TV. This property certainly has applications (for instance 
to celestial mechanics) which have not been explored yet. 

2.2 Paths of canonical transformations 

A striking result is that any continuously differentiable path of canonical 
transformations passing through the identity can be viewed as the Hamil¬ 
tonian flow determined by some (usually time-dependent) Hamiltonian. We 
will prove this in Proposition [2J and draw interesting consequences about 
the group of Hamiltonian symplectomorphisms. 

2.2.1 Symplectic and Hamiltonian isotopies 

We will call a canonical transformation / such that / = f f H for some Hamil¬ 
tonian function H and time t = a a Hamiltonian symplectomorphism. The 
choice of time t = a in this definition is of course arbitrary, and can be re¬ 
placed with any other value different from zero noting that we have / = f^ a 
where H a (z,t) = aH{z,at)\ the usual choice is a = 1. 

Definition 1 A symplectic isotopy is a one-parameter family (ft)t of sym¬ 
plectomorphisms depending in a C 1 fashion on t € M and such that fo = ha¬ 
lf each ft is a Hamiltonian symplectomorphism, then (ft)t Is called a Hamil¬ 
tonian isotopy. 

It turns out that each symplectic isotopy is a Hamiltonian isotopy, in 
fact the flow determined by some time-dependent H: 

Proposition 2 Let (ft)t be a symplectic isotopy. (i) There exists a Hamil¬ 
tonian function H = H(z, t ) such that ft = f^■ (ii) More precisely, we have 
( ft)t = (f t H )t with 

H(z,t) = — f a (j^ft o ft x (\z),z) dX. (7) 

Jo 

Equivalently: 

H{z, t) = — [ a (X H (f i ~ 1 (Xz),z))d\. (8) 

Jo 
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Proof. (Cf. Wang’s [73]; for a more conceptual approach using differential 
geometry see Banyaga [3]-) The ft being Hamiltonian symplectomorphisms 
we have Dft(z) T JDft(z) = J for every z € M 2n . Differentiating both sides 
of this equality with respect to t we get, omitting the variable z and writing 
ft = dft/dt , 

(Df t ) T JDf t + DtfjD(f t ) = 0 

and hence 

J{Df t ){Df t )~ l = [j(Df t )(D(f t )- 1 )] T . 

Using the chain rule together with the identity D(ft )~ 1 = D(ffl l ) we have 


D(Jf t of-') = DJ{f t of^) 


l T 


which shows that the Jacobian of J(ft ° / f _1 ) is symmetric. It follows from 
Poincare’s lemma that there exists for each f 6 / a function Hf such that 
J(ft o ffl 1 ) = d z Hf and one verifies that this function is explicitly given by 
the formula 

H t (z)= [ Jd z (f t o ff^Xz) • zdX. 

Jo 

Setting H(z,t) = — H t (z ) this is precisely formula ([7]). ■ 


2.2.2 Linear and affine flows 

When the Hamiltonian flow is linear, Proposition [2] yields an explicit for¬ 
mula: 

Corollary 3 Let ( St)t be a symplectic isotopy in Sp(n). There exists a 
quadratic Hamiltonian function H = H(z,t) such that St is the phase flow 
determined by the Hamilton equations z = Jd z H. The Hamiltonian function 
is the quadratic form 

H = ~JS t Sfl 1 z-z (9) 

where S t = dS t /dt. In particular, if S t = e tx with X € sp(n) (the symplectic 
Lie algebra) then H = — ^ JXz ■ z. 

Proof. We have ft o ffl 1 = StSfl 1 ] applying formula ([7|) we get 

H(z,t) = -J a(s t Sfl 1 {\z),z S ) dX ( 10 ) 
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which is ([9]), taking into account the linearity of a and St- ■ 
Writing St and its inverse in block matrix form 


S t = 





(Dj -Bj\ 
{-Cl Af ) 


( 11 ) 


(the second formula following from the identity StJSf = J ) we have 


JStS - 1 


(C t Dj-DtCT D t Aj — C t Bj\ 
\B t Cj - A f Dj A t Bj - B t Aj) ’ 


it follows from ([9]) that we have the explicit expression 
H = \{D t Cj - C t Dj)x 2 + (C t Bj - D t A[)px + \{B t A T t - A t Bj)p 2 (12) 


for the Hamiltonian function. 


Example 4 Consider the one-parameter family of matrices 

cost v(t) sinu;(t) 

— sinw(t) coscj(t) 

where uj is a C 2 function of time such that ip(0) =0. We have 



JStSf 1 


f-u(t) 0 \ 
V 0 ~c(t)J 


hence ( St)t is the flow of the time-dependent harmonic oscillator Hamilto¬ 
nian 

H{z,t ) = ^y-{p 2 + x 2 ). 

The case of affine symplectic isotopies is a straightforward extension of 
the result above: 


Corollary 5 Let ( St)t be a symplectic isotopy in Sp(n) and t \— > zt a C 1 
path in M 2n with zo = 0. The symplectic isotopy ( ft)t defined by ft = StT(zt) 
where T(zt) is the translation z i—» z+zt is the Hamiltonian flow determined 
by 

H(z,t ) = -^JS t Sfl l z ■ z + a (z, S t zt) - (13) 


10 



Proof. We have ft(z) = S z (z + zt) hence 


ft(z ) = Sflz + z t ) + S t z t = T{S t z t )S t T(zt)z 


hence ft = T(StZt)StT(zt). The inverse of ft being given by 


ft 1 = (StT(zt))- 1 = T{-zt)Si 1 


we thus have ftf t 1 = T(StZt)StS t 1 . Hence, by formula l]7|). 


H(z,t) 


j\(r(S t z t )S t Sfl\Xz),z) dX 
j a (^S t S^ 1 (\z), z^j dX — J a(S t z t ,z)dX 


which is (1131) in view of formulas ([9]) and (HOD , and taking into account the 
antisymmetry of the symplectic form. ■ 

The argument above can be easily reversed, yielding an explicit solution 
of the Hamilton equations for a Hamiltonian function of the type 


H(z,t) 


1 


M(t)z 2 + m(t)z. 


(14) 


Let in fact (St)t be the flow determined by the homogeneous part Hq(z, t ) = 
i M(t)z 2 ; we have (Hamilton’s equations) St = JM(t)St hence Ho(z,t) = 
— ^ JStSfl 1 z. Setting z t = Sfl l Jm(t), we can thus rewrite H in the form 
m■ Summarizing, the flow ( f^)t determined by (I14p is given by 


f t H = S t T(z t ) , zt= f S^Jm{t')dt'. (15) 

Jo 

Let us illustrate this by a classical example: the time-dependent driven 
harmonic oscillator. 


Example 6 Consider the Hamiltonian function 

H{z, t ) = ^(p 2 + x 2 ) + f{t)x. 

In view of Example [7] the flow ( S t )t of the homogeneous part is given by 


_ ( cosf l(t) sinfl(t)\ 

5 * sinfi(t) cos £l(t)J 

f t(t) = f u(t')dt' + fl(0). 

Jo 
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Setting m(t) = (/(f), 0).we have 

z t = (t') sin Q(t')dt',-f(t') cos Q(t)dt'^J 

and it suffices to apply formula m 


2.2.3 The flow determined by Hq + H\ 

Assume we are given a “master Hamiltonian” Hq, which we perturb by 
another Hamiltonian H\ (an archetypical example in the case n = 1 would 
be the choice Hq = p 2 / 2 and H\ = V ( x , t)). The following result shows how 
to calculate the flow (f f H )t of H = Hq + H\ knowing (f f H °) and (f t Hl )- 


Proposition 7 Let H = Hq + H\. We have 

ft 1 = ft*°ft* 1 with H{(z,t) = i7i(/f 0 (z),t). (16) 

Proof. Since ff 1 and fff° f^ 1 are both the identity on M 2n it suffices to 

show that the f-derivatives of f t H and f t °f t 1 are equal. We have, using the 
product and chain rules, 


d 




j t f t H °(ft H hz)) + Df?°{f?'{z))j t f?kz) 

X Ho U^ffhz)) + Df?\f?\z))X H{ {f?\z)). 


By definition, H\ = H\ o f f H ° hence, using ({3]), 


Df t H °(f t H hz))x Ht ( f r(z)) = Dfr(fr(z))x H{ ((fn- L frfr(z)) 




Hq t . 


eH 0 \-l 


pH 0 . 


= XH.UfUthV)- 


Summarizing, 

j t (ft Ho (ft H kz))) = X Ho (f?°f t H '(z)) + X Hl (f t H °f t H '(z)) 
= X Ho+Hl (f t H %f t H '(z ))) 


which we set out to prove. ■ 
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Assume, in particular, that H is separable, that is H = Hq + H\ with 
Hq{z) = T(p) and H\{z) = V(x). The flow (f f H °) is given by 

) = (x + tdpT (p ), p) (17) 

and we thus have 

H{(x,p) = V(x + td p T(p)) (18) 

so that (if 1 ) is obtained by solving the Hamilton equations 

x = td x V(x + tdpT (p)) DpT (p) (19) 

p = -d x V(x + td p T(p)) (20) 

where D 2 T(p) is the Hessian matrix (= matrix of second derivatives) of 
T calculated at p. When T(jp) = i p 2 is the ordinary kinetic energy these 
formulas reduce to the simple system 


x = td x V(x + tp) , p = — d x V(x + tp). (21) 


Setting u = x + tp this system is equivalent to the equations u + d u V(u ) = 0 
and p = u, that is to Hamilton’s equations for H. There are potential 
applications of this result to the theory of symplectic integrators (see e.g. 
Chorin et al. [T2], McLachlan and Atela [62] . Wang [73]). 

Let us illustrate these formulas on an elementary example. 


Example 8 Assume n = 1 and let Hq(x,p) = \p 2 and Hi(x,p) = | x 2 . 
The Hamiltonian flow determined by Hq is identified with the one- 

parameter group of symplectic matrices St = ^ ■ We thus have H\(x,p) 

tj(x' + pt) 2 and the Hamilton equations are x = (x + pt)t and p = — (x + pt) 
hence the flow determined by H\ consists of the symplectic matrices 


„H 


ft 1 = 


cos t + t sin t sin t — t cos t 


— sin t 


cos t 


( 22 ) 


and we thus have 


which is the flow of the harmonic oscillator Hamiltonian H{z) = ^(p 2 + x 2 ). 
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2.3 The group Ham(n) and its universal covering 

2.3.1 Products of Hamiltonian isotopies and Ham(n) 

Let (fH)t and (// v )t be Hamiltonian flows; we assume that they exist for 
t € It = [—' T, T\. It follows from Proposition [7] above that 

with H#K(z,t) = H(z,t) + K((f t H )~ 1 (z),t). (24) 

(/i^)" 1 = ft with H(z,t) = -H(f t H (z),t). (25) 

In fact, setting Ho(z,t) = H(z,t ) and Hi = K((f^ { )~ 1 (z),t) we have Ho + 
H i = Hj^K hence, applying formula (fTUl) . 

fH#K _ _ f H 0 fH{ _ fH rK 

Jt ~ Jt ~ Jt Jt ~ Jt Jt 

which is (1241) . Formula (1251) immediately follows noting that ft 1 = 
(for an alternative proof see Hofer and Zehnder [55j, p.144). 

We next assume that all Hamiltonians H are constant outside a com¬ 
pact subset /C of M 2n . We will call such Hamiltonian functions compactly 
supported (this is a slight, but innocuous, abuse of terminology). The vector 
field Xu determined by such a Hamiltonian function H is complete (see for 
instance Abraham and Marsden EU) and hence the flow (f^)t exists for all 
times t. Hamiltonian flows associated with such Hamiltonian functions are 
the identity outside the compact set /C. Hamiltonian symplectomorphisms 
coming from compactly supported Hamiltonians form a subgroup Ham(n) 
of the group Symp(n) of all symplectomorphisms; it is in fact a normal 
subgroup of Symp(n) as follows from the conjugation formula ([6]). That 
Ham(n) is a group follows from the two formulas (1241) and (125[) above. It 
turns out that Harn(n) is a connected group. To prove this it suffices to 
show that every / G Ham(n) can be joined to the identity /j by a path in 
Harn(n). But by definition of Ham(n) there exists a Hamiltonian function H 
such that / = fi\ the path we are looking for is just ( ft i )o<t<i■ It follows, 
using Proposition El that Ham(n) is the connected component of the group 
Symp c (n) of all compactly supported canonical transformations, i.e. equal 
to the identity outside a compact subset of M 2n (for detail see Hofer and 
Zehnder [55]).. 

2.3.2 The universal covering and a homotopy result 

Consider the universal covering group Ham(n) |55] of Ham(n); the elements 
of Ham(n) consist of homotopy classes (with fixed endpoints) of Hamiltonian 
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( 26 ) 


isotopies {f^)tei T and the group law is given by 

[/«"][/«*] = [/,"**] 

where [f t H ] is the homotopy class (with fixed endpoints) of the isotopy 
(//Oteir- The following result shows that the product in Ham(n) can be 
defined using either the relation (1261) , or by concatenation of Hamiltonian 
isotopies: 


Proposition 9 Let (f^)t and (f R )t be two Hamiltonian isotopies. (i) The 
paths {fF #K )o<t<2t 0 and {f R<>R )i<t< 2 t 0 where 


f ft K for 0<t<t 0 
1 ft-tjtf for t 0 <t< 21 0 


(27) 


are homotopic with fixed endpoints fi\ and f R ^ R = f R(> K ■ (ii) The group 
law of Harn(n) can thus be defined by 


[/f][/f] = \f t mK ]- 


(28) 


Proof, (i) Rescaling time if necessary, it suffices to consider the case to = \- 
Let us construct explicitly a homotopy taking the path {ft H fi K )o<t<\ to the 
path (,/f 0*) 0<t<i» that is, a continuous mapping 


h : [0,1] x [0,1] —» Ham(n) 

such that h(t, 0) = f R f R and h(t, 1) = f R<:>K for 0 < t < 1. Since we 
want to preserve endpoints during the deformation, we require in addition 
that h(0,s) = Id and h(l,s) = fi for all s € [0,1]. Define h by h(t,s) = 
aft, s)b(t, s ) where a and b are the functions 


aft, s ) 


I d for 


/ 


H 

(2t—s)/(2—s) 


0 < t < f 
for | < t < 1 


and 


We have aft, 0) 


b(t,s) = | 
= ft 1 , b(t, 0) 


/»/( 2 - 8 ) for 0 < i < 1 - i 

f( x for | < t < 1. 

= f t K hence lift, 0) = f t H f t K -, 


similarly 


hft, 1) 


fff for 0 < t < \ 
f-2t-lfl for \ < * < 1 
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that is h(t, 1) = f^ K . The relations h(0,s) = Id and h(l,s) = for 
0 < s < 1 are easily verified, ii) That the group law of Ham(n) can be 
defined by (051) follows from (i) and (f26l) . ■ 

One checks that the path (f t H( ' }K )o<t<t 0 is an isotopy corresponding to 
the Hamiltonian 


H0K(z,t) 


K(z , t ) for 0 < t < to 
H(z, t — to) for to < t < 2to 


(29) 


which is discontinuous at t = to- In fact, using formula ([7]) in Proposition [2] 
we have 

H0K(z,t) = - ^ a Hf t K o (f t K )-\\z),z) dX = K(z,t ) 

■Jo 

for 0 < t < to, and similarly 

H0K(z,t ) = - [ a (|(/£t 0 /£) o (/^ 0 /t^) _1 ( Xz),z) dX 
Jo 

= - / 1 ^(l/."o“(/" t »r 1 ( A-),^)dA 

Jo 

= H(z,t- t 0 ) 


for to <t < 2to- 


3 The Groups Mp(n) and IMp(n) 

The symplectic group Sp(n) has covering groups Sp g (n) of all orders q = 
2, ...,oo. Among all these the double covering Sp 2 (n) plays a very special 
role in mathematics and physics, because it can be faithfully represented 
by a group of unitary operators acting on L 2 (M ri ). This group is called the 
metaplectic group and denoted by Mp(n). We will loosely speak, as is usual 
in the literature, about the “metaplectic representation (i of the symplectic 
group”; one should keep in mind that, strictly speaking, n is a representation 
Sp 2 (n) —> Mp(n) of the double cover of Sp(n). 

3.1 Definition and main properties 

The symplectic group Sp(n) is a connected Lie group contractible to the 
unitary group U(n), hence 7Ti[Sp(n)] — (Z,+). It follows that Sp(n) has 
coverings n ? : Sp q (n) —> Sp(n) of all orders q = 2,..., oo. An essential fact 
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is that the double covering Sp 2 (n) has a faithful representation as a group 
of unitary operators acting on L 2 (M n ); this group is the metaplectic group 
Mp(n). 


3.1.1 Quadratic Fourier transforms and generating functions 

The group Mp(ro) is generated by the quadratic Fourier transforms S\v,m 
defined as follows: let 


W ( x , x') = ^ Px 2 — Lx ■ x + \Qx' 2 (30) 

be a real quadratic form where P = P T , Q = Q T , detL / 0 and set 
A m (lF) = i m ^J | det L\ where the integer m corresponds to a choice of 
argdetL. For if € <S(M n ) 

Sw,m^(x) = {2^m) n/2 A m (VF) [ e^ w ^ x ’ x '^{x')dx' (31) 

Jw 1 

(with the convention argi = 7r/2). In fact: 


Proposition 10 (i) Every S € Mp(n) can be written (non-uniquely) as the 
product Sw,mSw',m' of two quadratic Fourier transforms, (ii) We have 


Sw,mSw',m' = Sw",m" 


if and only i/det(P' + Q) / 0 and in this case 

P" = P- L t {P' + Q)~ 1 L 
L" = L'(P' + Q)- l L 

Q" = Q' -L'(P' + q)-\l') t ; 


(32) 

(33) 

(34) 


and the Maslov index of S\y",m" is given by 

m" = m + rri — Inert(P / + Q) mod 4. (35) 

(Hi) We have (SW,™) -1 = 5 'w',m' where W'(x,x') = — W(x',x) and m' = 
n — m, mod 4. 


Proof. For (i) and formula (1351) in (iii), see Leray [58], de Gosson 1231 [28] : 
formulas (f32l) ((Ml) are obtained by matrix multiplication using (f37l) . ■ 


Remark 11 Formula k35\) identifies the integer m with the Maslov index 
modulo 4 on Mp(n); see Leray Souriau [7 DJ , de Gosson 
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The covering projection II : Mp(n) —Sp(n) is defined by its action on 
the generators S\y,m'- 


Proposition 12 We have n(5'yp m ) = Sw where S\y G Sp(n) is generated 
by the quadratic form W, that is 


(x,p) 


S w (x',p') 


p = d x W{x,x') 
p' = — d x 'W ( x , x') 


(36) 


When W is given by (1301) we have the explicit formula 

Sw = (pL -1 Q - L t PL -1 ) ' ( 37 ) 

Equivalently, if 

£ , detP/0 (38) 

is a free symplectic matrix, then its generating function is 

W(x, x) = iDB-'x 2 - B~ l x ■ x + \B~ x Ax 2 . (39) 

For zo = (xo,po) G M 2n the Heisenberg-Weyl operator T(zq) is the time- 
one propagator of Schrodinger’s equation associated with the translation 
Hamiltonian z i —> (t(z,zq), that is, formally, T(zq) = e ~ ia ( z i z o)/fr where 



o(z,zq) = p ■ x 0 -po ■ x 


(40) 


where p = — ihd x and x is multiplication by x. It is hence a unitary operator 


on L 2 (R n ) whose action is explicitly given by 

T(z 0 )if(x) = e *(w-*-3W-*oty( x _ Xo ). (4i) 

The Heisenberg-Weyl operators satisfy the product relations 

T{z 0 )f( Zl ) = e^ Zo ^f( Zl )f(z 0 ) (42) 

T(z 0 + zi) = e-^ z °^f( Zo )f{ Zl ) (43) 

(the first formula follows from the second, interchanging zq and zi). Let 
S G Mp(n) and S = n(S'); we have the important intertwining formula 

ST (zo)S- 1 = T(Sz 0 ) (44) 
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which can be proven [28, [33] by checking it on the quadratic Fourier trans¬ 
forms S\v,m. generating Mp(n). We note that this formula is (incorrectly) 
taken in some texts as a definition of metaplectic operators; the irreducibil- 
ity of the Schrodinger representation and Stone-von Neumann’s theorem 
are invoked to motivate this “definition”. However, (1441) only defines a pro¬ 
jective representation of Sp(n), and not its double covering group; if one 
wants a true covering group one has to carefully determine the cocycle as¬ 
sociated with this projective representation (see Reiter [65] for a detailed 
analysis). In [23] we have given an analytical construction of this cocycle, 
and in [30] [37] [38] we use a topological and cohomological method, using 
techniques from symplectic geometry (the Leray index, and the notion of 
signature of a triple of Lagrangian planes, further developed in [28] [33]h 

3.1.2 The inhomogeneous metaplectic group 

Consider now the Heisenberg group H(n): it is M 2n x M equipped with the 
group law 

(z, t)(z ', t') = (z + z',t + t' + \a(z, z ')); 

the mapping p : H(n) —> U(L 2 (W 1 )) (the unitary operators on L 2 (R n )) 
defined by 

p(z 0 ,t 0 )ip = eTi to T (zo)V’ (45) 

is a unitary and irreducible representation of H n called the Schrodinger 
representation of EI(n). Defining the action of the symplectic group Sp(n) 
on HI(n) by S(z,t ) = ( Sz,t ) the group WSp(n) is the semi-direct product of 
Sp(n) and H n : 

WSp(n) = Sp(n) x H(n) 

with group law given by 

(S,u){S',u') = (SS',u(Su')) 

for u = ( z,t ) and v! = ( z',t ') (we mention that Burdet et al. [11] give a 
detailed study of WSp(n) and its generating functions). 

Let p : Sp 2 (n) —> Mp(n) be the metaplectic representation of the double 
cover of the symplectic group and let us define 

X : Sp 2 (?r) x H(n) —> U(L 2 (W n )) 

by X = P ® P- 

x{S,u) = p{S)p(u) , S G Sp 2 (n) , u = ( z,t ). 

We have: 
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Proposition 13 x unique unitary representation of the universal 

cover o/WSp(n) whose restriction to the Heisenberg group H(n) is p. 

We refer to Folland [2TJ (p. 196) for a detailed elementary proof of this 
result. The group defined by this extended representation is generated by 
the metaplectic operators and Heisenberg-Weyl operators, and is called the 
inhomogeneous metaplectic group IMp(n). It follows from the intertwining 
relation (144[i that every element of IMp(n) can be written T(zq)S (resp. 
ST(zo)) for some S G Mp(rt) and zq G M 2n . See Binz and Pods [5] for 
a detailed study of the Heisenberg and metaplectic representations with a 
detailed study of applications to optics, quantization, and held quantization. 

3.2 The path lifting property for Mp(n) 

Here is a fundamental property of the metaplectic representation, which is 
often undeservedly ignored in the physical literature. 

3.2.1 Lifting paths to the covering group 

To understand this, let us first state the path lifting property for covering 
spaces in its full generality (see e.g. Spanier ED). Let 7 r : Y —> X be 
a covering projection (Y and hence X are assumed to be connected) and 
7 : [a, b] —> X a continuous path such that 7 (to) = xq = 7r(yo) for some 
to G [a, b]. Then there exists a unique continuous path 7 : [a, b] — > Y such 
that 7707 = 7 and 7 (to) = Vo- Let us now choose X = Sp(n), Y = Mp(n) and 
assume that 0 G [ 0 , 6 ]. We choose for 7 a symplectic isotopy ( St)t in Sp(n) 
(thus S 0 = I). It follows from the path lifting property that there exists 
a unique path 7 of metaplectic operators St G Mp(n) passing through the 
identity in Mp(n) at time to = 0 and such that n(St) = St for every t G [a, b\. 
One shows [28l :33). using the fact that the Lie algebras of Sp(n) and Mp(n) 
are isomorphic, that the path 1 1 —> 7 (t) = St is C 1 if t 1 —> 7 (t) = St is, and 
that it satisfies Schrodinger’s equation 

i}Ss t = HS t , S 0 = 7 d (46) 

at 

where H is the quantization of the (time-dependent) Hamiltonian function H 
determined from (St)t using formula Q in Corollary [3j By “quantization”, 
we mean here the operator obtained by applying the Weyl correspondence 
to H (this will be detailed in Section [4.ID . Summarizing: 
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Proposition 14 Let H be a Hamiltonian function of the type 

H{z,t ) = ^M(t)z 2 

where M(t ) £ Sym(2n,M) depends in a C J (j > 2) fashion on t £ R. The 
solution of the Schrodinger equation 

_ 

ih-^ = Hfi , v>(-,o) = ^o a 2 (i n ) 

is given by if(x,t ) = St'if o(x) where (St)t is the unique path of operators 
St £ Mp(n) such that 11(5*) = St where ( Sf)t is the flow determined by 
Hamilton’s equations for H. 


Example 15 Let us illustrate this procedure on the harmonic oscillator 
Hamiltonian H(z) = \{p 2 + x 2 ). The flow (ft I )t consists of the rotations 

f H _ f cost sin A 
t \^— sinf cos tj ' 


Using formulas $39f) and \31\) the metaplectic operators with projection f t H 
are, for t ^ 7 rZ, the quadratic Fourier operators 


S^ifix) = 


i \ 1/2 i m r 
27 xih) | sin t J~ c 


exp 


i (x 2 + x /2 ) cos t — 2xx' 
h 2 sin t 


if{x')dx'. 
(47) 


One shows (Souriau m de Gosson that if we choose the value 

m = — [t/7r] ([■] the integer part function) for the Maslov index, then (Sf 1 ) 
is the lift of (f^)t to Mp(n), and is hence the solution of Schrodinger’s 
equation m> with Hamiltonian operator ^(—H 2 d 2 + x 2 ). 


3.2.2 Translation along a path 

Recall that the Heisenberg-Weyl operator T(zo) is the time-one evolution 
operator of Schrodinger’s equation 



a(z 0 ,z)^. 


where z = (x, —ihd x )- The following result (Littlejohn (59]) is an extension 
to the case where zq depends explicitly on t: 
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Proposition 16 Let 1 1 —> zo(t), zo(0) = zq be a C 1 path in the phase space 
M 2r L The solution of Schrodinger’s equation 


d'tp 

ih-7j£ = a(z, zo(t))ij} , ip(-,0) =ifo 

(48) 

is given by 

ip(x,t) = e^ x( ' t) f(z 0 (t))'tlJo(x) 

(49) 

where the phase x is real and given by 


If 1 , 

x(t) = — 2 y o a ( z o(t ),zo(t )) dt ■ 

(50) 

Proof. The idea is to write the evolution operator for the Schrodinger equa¬ 
tion (USD as the limit of a “time-ordered product”. Let to = 0 ,t\, ...,£jv = t 
be a sequence of successive times with A t = sup y \tj +1 — tj\ and set zq = 


zo(0), z\ = z(t \),zn = z 0 (t). Using the product formula (H3l) for the 
Heisenberg-Weyl operators we have 


T(z N - z N -i) ■ ■ ■ T(z 2 - z 1 )f(z 1 - zo)f(zo) 

( . N -1 

= exp ( 2h ^ <T ^ Zk + 1 ~ Zk ' Zk ) 

\ n k =0 

and one finds that in the limit N —>• oo ( i.e. At —>• 0) this product converges 
to the operator 

Ut = exp ^ a o(z(t'),zo(t'))dt , 'j T(z(t))- 

let us check that if = Ufi/j o indeed is the solution of (1381) . Writing (H9l) as 

if(x, t ) = e*^ x,t \p(x — xq (t)) 
the phase $ being given by 

<4>(x,t) = -± [ a(z 0 (t'),z 0 (t'))dt' - ^p 0 {t)x 0 (t) + p 0 (t)x 
Jo 

differentiation with respect to t yields 

ih^ = [<j(z 0 (t'), z 0 (tf))ip(x - x(t)) - ihxo(t)d x if(x - x 0 (t))] e 
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A similar calculation, differentiating with respect to the x variables shows 
that the right hand side is precisely a(jz, z(t))ip(x,t). ■ 

Suppose that z(T) = z( 0) = 0 and set j(t) = z(t) for 0 < t < T; then 

x(T) = — - j pdx — xdp = — j pdx 

and hence 

V>(x,T) = exp ( -jrf 7 pdx ) ip(x, 0); 

the initial and final states are the same up to a phase factor; this the occur¬ 
rence of a geometric phase shift, which is a phase difference acquired over 
the course of a cycle (“Berry’s phase” E3). 

Propositions [14] and [16] will be combined in Section 14.2.31 to explicitly 
solve the Schrodinger equation associated with Hamiltonians of the type 

H(z, t) = i M(t)z 2 + z 0 (t) • 2 . 

As noted in Burdet et al. CD, §6, there are difficulties in applying directly 
the path lifting property, but we will derive the result very simply, using a 
factorization result for evolution propagators. 


3.3 Mp(n) and Feichtinger’s algebra 
3.3.1 The Feichtinger algebra So(M n ) 


The metaplectic group Mp(n), and its inhomogeneous extension IMp(n), 
consist of unitary operators on L 2 (M n ); it is quite usual to perform practical 
calculations using the Schwartz space <S(M n ) which is dense in L 2 (M n ), and 
it is has become quite usual in quantum theory to use the latter as a natural 
“reservoir” for wavepackets. This choice is however inconvenient, and this 
for two reasons. The first is of a mathematical nature: 5(M n ) is a Frechet 
space, defined by an infinite number of seminorms, and this can make the 
proof of continuity properties for operators quite complicated. The second 
inconvenience is physical: the elements of <S(R n ) are C°° functions, and this 
is very restrictive from a physical point of view since its excludes de facto 
many realistic wavepackets (sharp pulses, step functions, etc.) such that, 
for instance, 


if{x) 


1 — |x| if |x| < 1 
0 if |x| > 1 


There is however an alternative choice to <S(M n ) which has none of the short¬ 
comings above. It consists in using the Feichtinger algebra Sq( M n ) (and its 
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generalizations, the modulation spaces M|(R n )), which was introduced by 
Feichtinger [181 [T9l [20] hr order to deal with localization problems in func¬ 
tional analysis arising in signal theory and time-frequency analysis. While 
signal theorists define £o(M n ) using “modulation operators”, it is better, 
keeping quantum mechanics in mind, to recast its definition in terms of the 
cross-Wigner transform 

W{if,f)(z) = (^hiT [ er^ v ' y ^{x + \y)4>*{x-\y)dy. (51) 

JR" 

Definition 17 The Feichtinger algebra So(M n ) consists of all if £ SfW 1 ) 
such that W(if, f) £ L 1 (M 2n ) for every cf £ 5(R n ), f ^ 0. The number 

IMU,So = W{ifA)\\F = [ \W(if,4>){z)\dz (52) 

JR 2 " 

is called the Wigner norm of if relative to the window f. 

A striking, but not immediately obvious, fact is that all the norms (1521) 
are equivalent when cf ranges over S( R ra ), and define a Banach space struc¬ 
ture on the Feichtinger algebra. The “windows” used in the definition of 
So(R n ) can themselves be chosen in S'o(R ri ): 

Proposition 18 Let both if and (f be in L 2 (R n ). (i) IfW(if,cf ) £ L 1 (R 2n ) 
then both if and <f are in ^(R"); (ii) We have if £ So(R n ) if and only if 
W(if,cf) £ L x (R 2n ) for one (and hence every) if £ S^R™). 

It immediately follows from this characterization of So(R n ) that: 

Corollary 19 A function if £ L 2 (R n ) belongs to So(R n ) if and only if 
Wif £ L 1 (R 2n ). 

Proof. In view of the statement (i) in the Proposition above the condition 
Wif £ L 1 (R 2n ) implies that if £ So(R n ). If conversely if £ So(R n ) then 
Wif £ L 1 (R 2n ) in view of the statement (ii) in the same Proposition. ■ 
The elements of So(R n ) are continuous functions; in fact: 

S 0 (R n ) C C°(R n ) n L 1 (R n ) n F{L 1 { R n )). (53) 

It easily follows from Riemann-Lebesgue’s lemma that each if £ Sc^R”) is 
bounded and that we have limi^^^ if = 0. The Feichtinger algebra is a 


24 


Banach algebra, both for pointwise product and for convolution. In fact, 
ip £ L 1 (M n ) and ip' £ So(M n ). Then ip * ip' £ So(M. n ) and we have 

\\^ * ^\U,S 0 ^ I MU 1 Ill'll*,So (54) 

for every window cp £ S(M. n ). Thus, if ip £ L 1 (R n ) and ip' £ So(M n ) then 
ip * ip' £ So(M n ), so that 

L 1 (R n )*S 0 (R n ) C S 0 (R n ). 

That Sq(R n ) also is closed under pointwise multiplication follows, taking 
Fourier transforms and using the fact that ip £ S'o(M n ) if and only if Ftp £ 

S 0 (R n )- 

3.3.2 Application to Schrodinger’s equation 

The fact that F(So(R n )) = 5o(M”) is a particular case of the following 
essential property: 

Proposition 20 The Feichtinger algebra So(R n ) is closed under the action 
of the inhomogeneous metaplectic group IMp(n). 

Proof. We recall the two following properties of the cross-Wigner transform 
[28 .. 33]: for every S £ Mp(n) we have 

W(SiP,ScP)(z) = W(iP,cP)(S- 1 z) 

for S = II (S'); for every zo £ M 2n 

W(f(z 0 )ip,f(z 0 )4>)(z) = W(ip, <p)(z - z 0 ). 

Since IMp(n) is generated by the operators T{zq) and S £ Mp(n) it suffices 
to show that if ip £ So(K n ) then T{zo)ip £ So(M n ) and Sip £ So(M n ). Taking 
definition (l52|) into account we have 

\\T{zo)ip\\ 4 ,,s 0 = [ \W{T{zo)ip,(p){z)\dz 
J M 2 ™ 

= [ \W(f(zo)ip,f(-z 0 )<p)(z-z 0 )\dz 

J R 2 ™ 

= 11^1 It(— 2 O )0,5 o 
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hence T{zo)ip £ So(M n ) since T(—zo)(p £ S(M n ). Similarly, 



\W(S^)(z)\dz 

\W(il),S- l (t))(S- l z)\dz 


S-i<t>,So 


and S~ l (p £ 5(R n ) hence Sip £ So(M n ) since the Wigner norms || • H^Sq and 
|| • | |^_i 0 So are equivalent. ■ 

It turns out (but we do not prove it here, see [46!, 33]) that So(M n ) is 
the smallest Banach algebra containing the Schwartz class <S(M n ) which is 
closed under the action of the inhomogeneous group IMp(n). 

An immediate important application of Proposition [20] is: 


Corollary 21 Let H be a Hamiltonian function of the type 

H(z, t) = i M(t)z 2 

where M(t) £ Sym(2ra,R) is a C' jf (j > 2) function of t £ M. Let ipo £ 
So(W l ). Let ip be the solution of Schrodinger’s equation 

ih W = ^ = ^°- 

The function ip(-,t) belongs to 5o(K n ) for every t £ R. 

Proof. It immediately follows from Propositions [14] and [20] ■ 

This result will be extended to inhomogeneous quadratic Hamiltonians 
in Section 14.2.31 


3.4 The Weyl symbol of a metaplectic operator 

Metaplectic operators map 5(M n ) —> 5(M n ) and it therefore makes sense to 
speak about their Weyl symbol. (We will review in some detail the notion 
of Weyl symbol of an operator in Section 14.11 to which we refer for the 
definitions used here.) 


26 




3.4.1 The symplectic Cayley transform 

Let Sp 0 (?r) be the closed subset of Sp(n) defined by 

Spg(n) = {S £ Sp(n) : det(5 — I) 0}. 

We have proven in [27; [29j [38] that: 

Proposition 22 (i) The mapping : S i —Mg defined by 
<S>(S)=M s = ±J(S + I)(S-I )- 1 

is an injection o/Sp 0 (n) into the set Sym(2n,R) of real symmetric 2 n x 2 n 
matrices. (ii)The inverse of is given by 

<h- 1 (M) = (M-ij )- 1 (M+ij) 

and we have <h(5 _1 ) = —<h(5). 

We call the mapping the symplectic Cayley transform and Mg the 
Cayley matrix associated with S £ Sp(n). 


3.4.2 The operators R U (S) 

Let now v be an arbitrary real number and define, for S € Sp 0 (n), an 
operator R U (S ) by the formula 


R, 


•(S) = {^) n iW\det(S-I)\ [ T(Sz)T(—z)dz. (55) 

JR 2n 


Using the formula (H3|) it is easy to rewrite (1551) as 

1 \ n 

1 \ l 


Ru(S) = 


T(z)dz. 


(56) 


2itHJ y|det(5-/)| J 

We have proven in [27] (also [291 [38] 1 the following result: 

Proposition 23 (i) Let S\y, m £ Mp(n) be such that II h (S\y,m) = 5 'w £ 
Spo(n). We have 

Sw, m — Rm— Inert W X x (Sw) (57) 

where Inert W xx is the index of inertia of the quadratic form x i— > W ( x , x). 
(ii) Every S £ Mp(n) can be written as a product Sw,mSw',m' with Sw , Sw' £ 
Spo(n) and if S = n^S) £ Sp 0 (n) then 

S = R m +m '+i sign($(S)+$(S')) ) ( 58 ) 

where sign(M) is the signature of the symmetric matrix M. 
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Remark 24 Formula {57 p identifies, modulo 4, the integers = m —Inert W xx 
with the Conley-Zehnder index of a path joining the identity I to S in Sp(n); 
see de Gosson f3(K IfSf /. 


Proposition [23] has the following consequences: 


Corollary 25 (i) The twisted Weyl symbol (sw,m )(7 of Sw,m £ Mp(n) (S\,y € 
Sp 0 (n) y ) is given by 


(sw,m)a(z) 


Am— Inert W xx 

i _X_ 

V|det(S w - I)\ G2h 


HSw)z 2 


(59) 


(ii) Assuming in addition that det(«SV + -0 / 0 the Weyl symbol sw of Sw,m 
is given by 


s\v,m(z) = 2"/ 2 


Inert W xx 

\/| det(SV + I)| 




(60) 


Proof. Formula (1591) is an immediate consequence of (1561) and (1571) in view 
of (1671) . To prove formula (1601) it suffices to note that sw,m is the symplectic 
Fourier transform (l 68 |) of ( sw,m)a (because the latter is involutive) and to 
apply the well-known Fresnel formula giving the Fourier transform of a 
Gaussian (see [28] . §7.4). ■ 

The case of the inhomogeneous metaplectic group easily follows: 


Corollary 26 The twisted Weyl symbol ofT(zo)Sw,m IMp(n) (with Sw £ 
Sp 0 (n) y ) is given by 


a a (z) 


pm —Inert W xx 

V 7 ! det(5 — I)\ 


e2h 


[$(S)(z-zo) 2 -cr(z,zo)\ 


(61) 


Proof. Using successively (1561) and (1431) we get 
1 


T(zo)Sw,m = 


2 vr h) ^|det( 5 -/)| 
1 i v 


2ttHJ y/\det(S-I)\ 


[ e ^ {s)z2 T(z 0 )T(z)dz 

JR 2 " 

[ e^ s > 2 e^ Zo ^f(z 0 + z)dz- 

JR 2 ™ 


formula ([6TD follows by the change of variables zq + z \—>■ z. ■ 

The Weyl symbol of a metaplectic operator S such that det(5 — I) = 
0 cannot be determined using the methods above; one has to use direct 
methods. 
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Example 27 Assume St 


1 t 
0 1 


for t € R. We have det(5t — I) = 0 


hence the formulas above do not apply. However, using the kernel formula 


m a direct calculation using Fresnel integrals leads to the following partic¬ 


ularly simple expressions for the Weyl symbol st of St and of its symplectic 
Fourier transform ( st) a ■ 



(62) 


4 Quantum Propagators and Isotopies 

Schrodinger’s equation is the quantum analogue of Hamilton’s equations; 
symplectomorphisms are replaced with unitary operators obtained from the 
Hamiltonian function by a quantization procedure. The process can be 
reversed, by dequantizing these unitary operators. This is consistent with 
Mackey’s view [0()| following which quantum mechanics is a refinement of 
Hamiltonian mechanics; (also see Marsden [63], Emch [16]). Mackey added 
that dequantization is a more fundamental process than quantization. 

4.1 The problem of quantization 

The problem of how to “quantize” properly a “classical observable” is still 
an open one. Mathematically, it amounts to finding a pseudo-differential 
calculus suitable for applications to physical problems. 

4.1.1 Shubin and Weyl correspondence 

Let us review some basic concepts from pseudo-differential theory following 
Shubin [69]. Let A be a continuous operator <S(R n ) —> 5'(R n ); in view of 
Schwartz’s kernel theorem m there exists a distribution K € x R n ) 


such that 


f))). 


((•, ■) and ((•, •)) the distributional brackets on R n and R n xl", respectively). 
Writing formally 



Let now r be an arbitrary fixed real number; the r-symbol a T £ S'(R 2n ) of 
A is defined by the Fourier transform 
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One shows [281 |33l [69] that A can be (formally) written as a pseudo¬ 
differential operator 

Aif(x) = ( 2 ^ 1 )" / [ e^ p ' {x ~ y) a T {{ 1 - t)x + ry,p)ip(y)dydp. (64) 

J J R 2 " 

We will write A = Op r (a T ) and call A the r-operator with symbol a T . 
In particular, the choice r = \ corresponds to the Weyl operator A = 
Op™(a 1/2 ); Weyl symbol is the distribution a £ tS^R 2 " - ) given by 

a(x,p)=[ e~^ p ' v K{x + \y,x-\y)dy (65) 

it" 

(see e.g. Shubin [69j . de Gosson |2Bj 331)- We will write A = Op™(a); the 
operator A is then given by 

M{x) = {^nT [ e~& y a{\(x + y),p)if{y)dpdy] (66) 

I1 2 " 

one proves that A is self-adjoint if its symbol a is real (this property is not 
shared by the operators (1641) for t ^ 

One shows that the following harmonic representation of A holds: 

Op “(a) = ( 2 ^)" [ a a (z)T{z)dz (67) 

where a a is the symplectic Fourier transform of the symbol a; formally 

“AV = (if [ e-i°^a(z')Az'. (68) 

JM. 2n 

The distribution a a is called the twisted (or covariant) Weyl symbol of A. 

A characteristic property of Weyl operators is their symplectic covari¬ 
ance; it is the quantum analogue of the property ([6]) of Hamilton’s equations. 

Proposition 28 For every S £ Mp(n) with S = n(S') we have 

SOp w {a)S~ 1 = Op w (aoS~ 1 ). (69) 

Conversely, if A = Op T (a r ) is such that S Op T (a T )S~ l = Op r (a T oT 1 ) 
then we must have r = k; i.e. A is a Weyl operator. 

A similar property does not hold for arbitrary Shubin r-operators. See 
however de Gosson [35j for partial symplectic covariance results. 


30 


4.1.2 Moyal product and twisted convolution 

The Moyal (or star) product a *n b of a, b G <S(M 2n ) is defined by 

a * h b(z ) = f ’ zin> a(z + \z')b{z - \z")dz!dz"\ (70) 

jr 2 " 

this formula can be alternatively written 

a* h b(z) = ( 4^) 2 " / e- 2 i d < z ’ z '^a{z')b{z")dz'dz"- (71) 

JR 4 " 

where da is the coboundary of the symplectic form: 

da(z, z', z") = a(z, z') — a(z, z") + a{z ', z"). 

The twisted convolution a#6 i defined by 

a#b = F a (F a a *n F a b); (72) 

explicitly: 

a # b ( z ) = (i mT [ e^ z ’ z '^a(z - z')b{z')dz’ (73) 

JR 2 " 

The relation of these notions with Weyl operator theory is the following 

EH S3 EH: 

Proposition 29 Let A = Op™(a), B = Op™(6) and C = AB. We have 
C = Op™ (a *n b) and (a *n b) a = a#6. 

The papers by Hansen m and Littlejohn [59] contain nice reviews of 
the topic; also see the evergreen book by Folland m, and de Gosson |28l 155] 
where the closely related topic of deformation quantization is also discussed. 

4.1.3 Born—Jordan quantization 

It is the property of symplectic covariance (1691) . and the fact that A = 
Op™(a) is self-adjoint if a is real, which are at the origin of the privileged 
role played by Weyl operators, especially in quantum mechanics. The rub 
comes from the fact that it is not quite obvious that Weyl quantization is 
really the right choice in physics; for instance we have shown in [36] that if 
one wants the Schrodinger and Heisenberg pictures of quantum mechanics to 
be equivalent (which is assumed to be true by most physicists), then one has 
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to use the so-called Born-Jordan quantization scheme m we have studied 
in Easu; the latter is obtained by averaging the r-symbol (|63l) over the 
interval [0,1]: 


/' 


a B j{x,p) = / a r (x,p)dT. 


(74a) 


In the case of polynomials, the distinction between Weyl and Born- 
Jordan quantization is known since the early years of quantum mechanics; 
for instance the Weyl correspondence is the prescription (in dimension n = 
1 ) 


x q 


Weyl 1 
* 2 1 


£ 

£=0 


ff- e x r f/ 


(75) 


while Born-Jordan imposes the apparently less weighted rule 


x s q r 


BJ 


S + 1 


1=0 


- s ~ e x r f/ 


(76) 


It turns out, however, that Weyl quantization and Born-Jordan quantiza¬ 
tion coincide for large classes of operators; for instance for all physically 
interesting operators associated with Hamiltonians of the type T(p) + V{x) 
or, more generally 



(see de Gosson ES]), so one can safely conclude that within the framework 
of the present paper it really doesn’t matter whether one uses the Weyl or 
Born-Jordan quantization rules. 

4.2 Paths of unitary operators 

Hamilton’s equation govern the time evolution of classical systems; in non- 
relativistic quantum mechanics Schrodinger’s equation plays a similar role. 
While points in the classical phase space are propagated using canonical 
transformation, the evolution of quantum-mechanical wavefunctions is ob¬ 
tained using paths of unitary operators. 
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4.2.1 Schrodinger’s equation 

Let H be the Weyl quantization of a Hamiltonian function H; in physics H 
is typically of the type 


n 

H = 12 - 4/(M)) 2 + V(x,t) 


3 = 1 

where the potential functions Aj and V satisfy some smoothness conditions. 
When H is of the type above this operator is given by 


1 


d 


1=1 3 


in both the Weyl and Born-Jordan quantization schemes. The Schrodinger 
equation corresponding to the Hamiltonian function H is 

^ = t) , V’('T') = P (77) 

where the initial Cauchy datum ip' is a function (or distribution) belonging 
to some suitable function space (Yajima [751, 76] has introduced a class of 
distributions which guarantees the existence of the solutions). Assuming 
that the solution ip of (1771) exists and is unique for t in some interval I = 

[t' — T,t' + T] we can write ip = U^ t ,ip' where U^ t , (the evolution operator, or 
propagator) is a unitary operator on L 2 (R n ). The Chapman-Kolmogorov 
property 

U t H t ,U t V t „ = U t H t „ , U t H t = I d (78) 

holds whenever U^, Ujf t „ exists. When t' = 0 we write ip' = ipo, It = 

[—T, T], and U t H 0 = Up. Notice that the semigroup property U^Up = Up +t , 
only makes sense when H is time-independent. 

Since Up t , is a continuous operator 5(R ra ) —> L 2 ( R n ) C 5'(R n ), it is a 
bona fide Weyl operator; the determination of its Weyl symbol is however in 
most cases a cumbersome exercise (to say the least); Berezin and Shubin [7] 
give an expression for the Weyl symbol using a Feynman path-type integral, 
which is not very tractable in practice. 

Suppose now that the Hamiltonian function is a (time-dependent) quadratic 
form in the position and momentum variables; such a function can always 
be written as 

H(z,t) = \M(t)z 2 (79) 
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where t i— > M(t ) is a C ] (j > 2) mapping R —> Sym(2n, R). In that case 
we have U^ t , = S^t 1 where the two-parameter family is constructed 

as follows: setting = U t H {) we have U^ t , = U t H . Consider now the 
(time-dependent) Hamiltonian flow (f^)t', it consists of a C 1 one-parameter 
family of linear mappings S t £ Sp(ra) passing through the identity at time 
t = 0. Using the path lifting property previously studied, to (ft H )t = {St)t 
corresponds a unique family (S t ) t of metaplectic operators passing through 
the identity at time t = 0 hence we have 

U^ t i = £ Mp(n). (80) 


4.2.2 Quantum isotopies 

In Proposition[2]we showed that each symplectic isotopy is actually a Hamil¬ 
tonian flow. We are now going to prove a quantum analogue of this property. 
We begin by defining the notion of quantum isotopy, which is the operator 
analogue of a symplectic isotopy. 


Definition 30 A quantum isotopy is a C 1 one-parameter family {Ut)t of 
unitary operators on L 2 (M ri ) having the following properties, (i) JJq is the 
identity operator: Uq = I a; (ii) There exists a dense subspace D of L 2 (W n ) 
such that 



lim 

—^0 


Ut+At'f’ - Ujlf 

At 


exists for every £ D. We call D the domain of the quantum isotopy ( Ut)t ■ 
When Ut € Mp(n) for every t £ R we call ( Ut)t is a metaplectic isotopy. 


We will also use the following notation: for t, t' £ R we set U t u = UtU t , 1 
hence U t ,o = U t , U t ,t = Id, and 

Ut,t'Up r = U t ^ , (tvr 1 = £/*',* (81) 

(c/. the Chapman-Kolmogorov law (1781) 1. 

A quantum propagator (U^)t is de facto a quantum isotopy: we have 
U^ 1 = /a and the relation 


fuf ) * = -M 


1 

ih 


holds for if £ <S(R n ). Writing U t = 4fUt we have, more generally: 
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Proposition 31 Let ( Uf)t be a quantum isotopy with domain D. (i) The 
(time-dependent) operator 

H = ihUtUf 1 (82) 

with domain Dg = D is self-adjoint; (ii) the function f) = Ut'ifo is a solution 
of Schrodinger’s equation 


ih ^t = ^ = 

for every 'f/’o £ Dg; (in) When ( Ut)t is a metaplectic isotopy ( St)t, then H 
is explicitly given by 


1 


H = --JS t Sf l z-z 


(83) 


where ( St)t is the symplectic isotopy obtained by projecting ( St)t on Sp(n). 

Proof, (i) The operator Ut is defined on L 2 (M”) and the domain of H is 
hence that of Ut■ Let us set 


ih w i ih 


-1 


^At = -fe(Ut+At - Ut)U^ = —{U t+At ,0 - U t ,o)U t ~ fi 

that is, since (Uf q 1 )* = Ut, and using the rules (18TT) . 

ih , , ih , 

Aa t = U t+ At,o ~ Ut,o)Uo,t = -^(Ut+At,t ~ I)- 

It follows that for £ D, 


(84) 


lirn {A At if\<t>) = ih(UtU t l ip\<f>) = {Hif\4>). 

At-tO 


(85) 


Similarly, taking into account the equality Uf +At t = Ut.t+At, and using again 
the rules (1511) . the adjoint of A At is given by 

ifl X Jlj 

A*a t = -~^{U t ,t+At - I) = - U t+ At,o)Uo,t + At 


that is 


ih 


A* At = —(U t+ At-U t )U- + \ t . 


( 86 ) 


It follows that for £ D we have 


lim {if\A At <j)) = ih(4>\U t U t f) = (tp\H<j>). 

At—>0 
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Since = (V’I^a^) We ^ ave th us proven that (Hip\(j)) = (ip\Hf) 

for all ip, 4> € D, which shows the self-adjointness of the operator H. (ii) 
immediately follows from the definition (1821) of H. (iii) Let us set St = II (St)', 
the one-parameter family (St)t is a symplectic isotopy in Sp(ra), and is thus 
(Corollary [3]) the Hamiltonian flow determined by a quadratic Hamiltonian 
function, the latter being given by formula (j9j) : 

H(z,t ) = —JStS^z ■ z; 

the expression (f83l) of H since {St)t is just the lift to Mp(n) of the symplectic 
isotopy (S t )t- ■ 

Stone’s theorem allows us to somewhat weaken the assumption on ( Ut)t 
when it is in addition a one-parameter group: 

Corollary 32 Let ( Ut)t be a strongly continuous one-parameter group of 
unitary operators on L 2 (M. n ): lim t^t 0 Utif = Ut 0 ip for every if € L 2 (W 1 ), 
and UtUt' = U t +t' for all t,t' € R. Then ( Ut)t is a quantum isotopy, and we 
have Ut = e~ lh ^^. 

Proof. By Stone’s theorem [SI El EH there exists a unique self-adjoint 
operator H (the infinitesimal generator), densely defined in L 2 (M n ), whose 
domain contains <S(M”), and such that Ut = e ~ lhH ^This shows that the 
path 1 1 —> Ut is C 1 and hence a quantum isotopy since ihUtUjf 1 = H. ■ 

4.2.3 A factorization result 

We are going to prove the analogue of Proposition [7] about Hamiltonian 
systems when Hq is a quadratic polynomial (Weinstein [13]) proves a similar 
result in the particular case Hi = V(x,t)). 

Proposition 33 Suppose Hq is a quadratic Hamiltonian of the type m 
Then the propagator i s given by the formula 

jjHo+Ih = JjHojjmoSt = g tU H ioSt 

where ( S t )t = ( ft I °)t is the flow determined by the Hamilton equations for 
Hq (and thus St E Sp(n)J. 
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Proof. Let us set V) = U^° U ^' oSt . We have, using respectively the product 
rule and metaplectic covariance formula (|69l) for Weyl operators, 

= H 0 U t HloSt + S t HT°StU t HloSt 
= H 0 U t HloSt + StiS^H^U^ 

= H 0 U^ loSt + HiU? loSt 


which proves formula (1871) . ■ 

The assumption that Hq is quadratic is essential in the proof of formula 
(|87li . since the flow (f^°) then consists of symplectic matrices, allowing us 
to use the symplectic covariance formula H\ o St = Sfl 1 HiSt- There is no 
analogue in the case of non-linear flows (see de Gosson ISj). However, we 
conjecture that it is possible to obtain asymptotic equalities in the limit 
h —> 0 (with arbitrary accuracy) using the method developed by Lasser and 
her coworkers |22l [56], and which yield a refinement of Egorov’s theorem 
m on transformation properties of pseudo-differential operators. 


Example 34 Assume that Hq = ^p 2 and Hi = V(x,t). Set H = Hq + H\. 
Then U = SflJY* where St is the free particle propagator 


s,<p(x) 


( * 'VV f exp hh^)!l 

\2mh\t\J J-oo [K 


if(x')dx' 


( 88 ) 


with m = 0 for t > 0 and m = 1 for t < 0, and Vflzfl ) = V(x +pt,t) (the 
latter being obtained by applying the path lifting method to Hq) since the flow 
determined by Hq consists of the linear mappings St : (x,p) i —> (x + pt,p). 


As an illustration let us work out in detail the case of the harmonic 
oscillator H(x,p) = \{p 2 + x 2 ) considered in Examples 151 and [T5l It gives 
rise to a group (U f H )t of unitary operators given, for if € <S(R) and t ^ 7rZ, 
by the formula fljTp . Let us set Hq(x,p) = \p 2 and Hi(x,p) = \x 2 . The 
evolution operator U t H ° is the quantization of the one-parameter group of 

St = ^ ^ and thus consists of the metaplectic operators defined, for 

if € <S(M) and t / 0, by (j88]i . We have H\{x,p) = ^(x + pt ) 2 , and the 
evolution operator JJ^ loSt is thus the quantization of the linear flow (|22[) . 
given by 



cos t + t sin t sin t — t cos t 
— sin t cos t 
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-H\ 


the symplectic matrices f t 1 are free for t p 0 and their generating functions 
(I3U1) are 


W (x, x', t) = — 


1 


sin t — t cos t 


1 2 , / , 1 / 2 / , , , ■ 

—x cos t — xx + —x (cosr + tsmrj 


It follows, in view of formulae (I3T1) and (1391) . that jjp loSt cons i s t s Q f the 
metaplectic operators defined by 

1 \ 1/2 




Hi_o St 


( 1 \ jm roo 

—— ) —/ e^ w ^ x,x ^(x')dx' (89) 

2irih J sint — t cost| J-oo 

for an adequate choice of the Maslov index m. Now U^°U^ loSt is the 
product of two metaplectic operators of the type ([3T]) . namely Sw,mSw',m' 
with P = L = Q = t -1 and 

cos t 1 


P' = 


sin t — t cos t 


,L’ = - 


The sum 


P' + Q = 


sin t — t cos t 
sin t 


, cos t + t sin t 
sin t — t cos t 


sin t — t cos t 

only vanishes for t £ nZ, hence formulas (f32]h (j33j), and (l34]l in Proposition 
E2 yield, for t £ 7rZ, U f H(> U f Hl ° St = S\y", m " with 

P" = Q" = ^1 ; L " = _L. 

sin t sin t 

We thus have U?°U? loSt = UP. 


4.2.4 Inhomogeneous quadratic Hamiltonians 

The results above allow us to prove the following extension of Propositions 
[20] and [TJ] to inhomogeneous quadratic Hamiltonians, thus solving a problem 
unsuccessfully addressed in Burdet et al. CD. §6. We first remark that every 
in every such Hamiltonian 

H(z, t) = -M ( t)z 2 + m(t) ■ z 

with M(t) € Sym(2n,M) and Zo(t) £ M 2ri , both being C(j > 2) functions 
oftGR, the inhomogeneous term can be written as a Hamiltonian of the 
type a(z,z(t )) where 

z(t) = J f m(t')dt'+ z( 0). 

Jo 
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Proposition 35 Let H be a Hamiltonian function of the type 

H(z, t) = -M(t)z 2 + a(z, z(t)). (90) 

(i) The solution of the corresponding Schrodinger equation 

d ib '' 

ih ~dt = H ^ ’ ^' ,0 ) = 
is given by if = Uj^ifo with 

U^ S t T(u(t) , u(t) = f Sfi lz (t')dt' (91) 

Jo 

where (St) is the evolution operator for the Weyl quantization of Ho(z,t) = 
\M(t)z ■ z and the phase x is given by the formula 

X ^ = ~\J a i s t'z(t'),z(t'))dt'. (92) 

(ii) If if o G So(K n ) (the Feichtinger algebra), then if (-,t) G So(M n ) for every 
t £R. 

Proof, (i) Write H = Ho + H\ where Ho(z,t) = | M(t)z ■ z and H±(z,t) = 
a(z,z(t)). In view of formula (1571) in Proposition 1551 we have _ 

StU^ 10 t where St is the evolution operator for Schrodinger’s equation with 
Hamiltonian operator Hq and St = n(S)) (Proposition 1141) . On the other 
hand, we have 

(Hi ° S t )(z, t ) = a(S t z , z(t)) = a(z, S^zft)) 

hence 

Hi o S t = a(z, S^zft)) = a(z, u(t )) 

where we have set u(t) = f ( - Sf 1 z(t')dt'. It follows from Proposition [T6l 
formula (1751) . that we have 

U t Hl ° St = exp cr(z(t'),u(t))dt'^j T(u(t)) 

= exp z(t))dt^j T(u(t )) 
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and hence 


U t H = exp J (r(S t 'z(t'),z(t))dt , '\ S t T[u{t )) 

which is (1911) . (ii) Follows from formula (1911) using Proposition [20] and ob¬ 
serving that, up to the unimodular factor e'ALA the propagator is in 
the inhomogeneous metaplectic group IMp(n). ■ 

Part (ii) of the Proposition above can be interpreted by saying that 
the “phase space concentration” of an initial wavepacket is preserved when 
the quantum propagator arises from a (homogeneous or inhomogeneous) 
quadratic Hamiltonian; we will discuss an extension of this result in next 
Section. 

4.3 An extension of the metaplectic representation 

In the 1980s Weinstein m suggested a way to extend the metaplectic rep¬ 
resentation by constructing solutions to Schrodinger’s equations with non¬ 
quadratic potentials. Weinstein’s “generalized metaplectic operators” have 
recently been rediscovered in m where the authors define and study an 
algebra of Fourier integral operators and their action on modulation spaces. 
They assume that the Hamiltonian is a quadratic function perturbed by 
a symbol belonging to certain classes of modulation spaces, the so-called 
Sjostrand classes |67’ ) 68]. 

4.3.1 The Sjostrand class 

The next application is a regularity result extending (ii) in Proposition [35] 
hereabove. Let us first introduce a convenient class of symbols (Sjostrand 
(67, 68], Grochenig [761135] ): 

Definition 36 The Sjostrand class M°° ; 1 (M. 2n ) is the complex vector space 
consisting of all a G S'(M? n ) such that 

sup \W(a, b)(z, -)| G L 1 ^ 2 ") 
zeR 2n 

for all b G 5(M 2n ); here (z,?) i —> W(a,b)(z,<j) is the cross-Wigner trans¬ 
form on M 2n x R 2n . 

The Sjostrand class is a Banach space for the topology defined by the 
equivalent norms || • \ \ b , M oo,i defined, for 4> G 5(M 2n ), by 

= / sup |W(a,4>)(z,C)K; 

J R 2n zeR 2n 
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the Schwartz space 5(M 2n ) is dense in M 00,1 (R 2n ). We also mention that 
M 00,1 (M 2n ) is invariant under any linear change of variables: if a € M 00,1 (M 2n ) 
and M G GL(2n, M), then a o M G M 00,1 (M 2ri ). 

In our context, the interest of this symbol space comes from the two 
following properties: 

Proposition 37 Let a G M 00,1 (R 2n ). (i) The Weyl operator A = Op w (a) 
is bounded on the Feichtinger algebra So(M. n ). (H) M 00 ’ 1 ^ 2 " - ) is a Banach 
algebra for the star-product *n: if b G then a b G M°° ,1 (M. 2n ), 

and for every $ G <S(R 2n ) there exists a constant C§ > 0 such that 

I l a *h b \< C$11a,| 1HI^.M 00 ’ 1 (93) 

for all a*fib e M°° ,1 (M 2ri ). 

We refer to Grochenig |36( @8] for a proof of these properties. 

4.3.2 A perturbation of the metaplectic representation 

In a recent paper Cordero et al. H33 establish the following result, which 
considerably extends Proposition 1201 

Proposition 38 Let a G M°°’ 1 (M 2n ) and set 

H(z,t) = ^M(t)z 2 + a(z,t). (94) 

Consider the associated Schrodinger equation 

= ^ = ^°' 

We have ip(-,t) G <So(R n ) for every t G R if and only ififo £ <So(M n ). 

This result is very interesting, because it shows that the perturbation of 
a quadratic Hamiltonian by an operator associated to a Sjostrand symbol 
does not affect the phase-space concentration of the solution to Schrodinger’s 
equation. Such perturbations can thus be viewed as “small” from the phase 
space perspective. Wether a similar result holds for larger classes of “per¬ 
turbations” is an open problem. 

It turns out that Proposition [38] can be obtained from the factorization 
result (1B71) in Proposition [33] We are going to sketch the argument here; the 
full details will be published elsewhere. Let us write H = Hq + H\ where 


41 


Hq = \M(t)z 2 and H\ = a; in view of ([571) we have = StUf oSt where 
(St)t is the metaplectic isotopy determined by Hq- In view of the closedness 
of So(M n ) under the action of Mp(n) it is thus sufficient to show that Uf oSt 
maps So(M n ) into itself. Since the Sjostrand class is invariant under any 
linear change of variables, the proof of Proposition [38] thus boils down to 
the simpler statement: 

Lemma 39 Let H € Af 00 » 1 (R 2n ). Then U t H : S 0 (R n ) —> S 0 (R n ) for every 
t € R. 

Sketch of the proof. For simplicity, we assume that H is time-independent; 
then one shows that the Weyl symbol of U f H can be written 

as a series 

u t = I + H + * h H + * n H * h Hf + • • • 

which is convergent in view of the inequality (1931) in Proposition [37] Hence 
ut £ M°° ,1 (R 2n ), and it then suffices to use part (i) of the same proposition. 

5 Discussion 

Phase space picture of quantum mechanics There is dissymmetry 
between Hamiltonian mechanics and quantum mechanics in its usual for¬ 
mulation (the Schrodinger representation, which we have been considering 
here). Hamilton’s equations are expressed in terms of a phase space, en¬ 
dowed with a symplectic structure, while Schrodinger’s equation leads to 
evolution operators defined on functions on configuration space. This dif¬ 
ficulty is actually not a real one, because one can replace the standard 
Schrodinger equation by the phase space equation 

= H (x + \ihd p ,p — \ihdx) 'P 

obtained from the Hamilton function H using Bopp quantization instead 
of Weyl quantization (see [sa sa su. This equation allows to define a 
propagator U f H acting on functions defined on phase space, which is re¬ 
lated to the usual propagator Ujf 1 by an intertwining relation. The phase 
space Schrodinger equation is actually a variant of deformation quantization; 
this approach might actually be, both mathematically and physically, more 
fruitful than the traditional Schrodinger approach, since it views quantum 
mechanics as a deformation of classical mechanics. 
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Epistemic considerations One cannot help being struck by the beauty 
and harmony of the structures underlying Hamiltonian flows and isotopies, 
and their operator counterparts, quantum isotopies. Both are governed by 
almost similar - one is tempted to say identical - algebraic laws correspond¬ 
ing to underlying physical processes that are well illustrated by Bohrn and 
Hiley’s “implicate order” [8]. As hinted at in the Introduction (also see 
the discussion in Section 13.2[) , one actually passes from Hamiltonian flows 
to their quantum counterparts by “lifting” the first to a covering group; 
since the most obvious occurrence of this property can be seen on the meta- 
plectic representation, we have called elsewhere 1261 SI] metatron the entity 
whose evolution is governed by the quantum motion induced by the quan¬ 
tum propagator. In that philosophical view, Hamiltonian flows can be seen 
as “unfolding” reality; Hiley’s work [53] might well help shed some light on 
the deeper algebraic structures; in Hiley’s language our product formulas 
could be interpreted as algebraic unfolding processes underlying a “holo- 
movement” (also see de Gosson and Hiley [40] for a musical comparison). 
There is however a fundamental difference between Hamiltonian and quan¬ 
tum mechanics, already visible at the level of mathematics: the first is local, 
while the second is non-local, global (c/. the Moyal product (|70l) where one 
integrates over the whole space at time t ). These issues of classical versus 
quantum processes will be further discussed in [ ill . 
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